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Abstract 



We calculate, numerically, the residence times (and their distribution) 
of a Brownian particle in a two- well system under the action of a periodic, 
saw-tooth type, external field. We define hysteresis in the system. The hys- 
teresis loop area is shown to be a good measure of synchronization of passages 
from one well to the other. We establish connection between this stochastic 

synchronization and stochastic resonance in the system. 
PACS numbers: 82.20.Mj, 05.40. +j, 75.60.Ej 
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Nature, presumably, exploits the phenomena of stochastic resonance (SR) to its advan- 
tage [|T],0 to tune in to a desired signal. However, its discovery, initially as a theoretical 
invention to explain recurrence of ice ages |||]], is barely a decade and half old. SR is a 
nonlinear phenomena wherein an input noise is partially rectified by the system to obtain 
enhanced output signal at the input signal frequency. It is reflected as a peak in the output 
signal-to-noise ratio (SNR) as a function of input noise strength. The nature of output signal 
in the nonlinear system, however, depends on the combined effect of the input signal and the 
input noise. Since the input signal, in general, is often considered to be weak (subcritical), 
in the absence of noise no output signal is obtained. The output signal considered is an 
averaged effect of a large number of realizations of fluctuating forces (noise) in combination 
with the input signal. The output signal, therefore, may not have the same form as the input 
signal. However, one expects the output signal to have its dominant time scales of variation 
nearly synchronous with those of the input signal. The degree of synchronization, however, 
will depend on the noise strength as well as the input signal parameters. It is, therefore, 
interesting to study how the output signal gets synchronized with the input signal, and, 
in particular, whether SR occurs exactly when there is maximum synchronization between 
the two signals. In this work we show, by calculating the residence time distributions in 
a periodically forced two-well potential, that there is a close connection between these two 
phenomena. 

We define hysteresis in a two-well system, as explained below, such that its loop area is a 
good measure of synchronization of passages from one well to the other under the influence of 
a Gaussian white noise f(t) when subjected to an external periodic field h(t). Just in order 
to make the explanation simpler and more transparent we define hysteresis loop from the 
distribution of first-passages from one well to the other. ( The explanation may subsequently 
be simply carried through to the main part of this work wherein residence time distribution 
is used.) The procedure is illustrated in Figures l(a-d). p(r) is the distribution of first- 
passage-times, r, from one well to the other (Fig. lb) obtained by solving the overdamped 
Langevin equation 
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numerically (for a large number of realisations), where 

U(m) = --m 2 + -m 4 - mh{t) (2) 

represents the two-well potential in the presence of an external field h(t) which is periodic in 
time t. The fluctuating forces f(t) satisfy: (/(£)) = and (f(t)f(t')) = 2D5(t — t'), where 
D is the strength of the noise. The sequence of r's is then used to calculate the distribution 
p(h(r)) (Fig. lc) of field values h{r) at which first-passages occur from the right well to the 
left one. From p(h(r)) we obtain the upper half of the hysteresis loop M(h) (Fig. Id), 



M(h) . L /•/"■ 



• 2 / p(h')dh', (3) 
h c Jh 

the other half being obtained by symmetry ||. Here h is the amplitude and \h c \ is the 
minimum value of the field h(t) (Fig. la) beyond which one of the two wells of U(m) 
disappears. 

From Fig. lb we see that the peaks of p(r) occur periodically and centered each time 
around which h(t) = — h at which the potential barrier of passage from the right to the left 
well is the least. Now, perfect synchronization of passages would mean sharp (^-function) 
periodic peaks in p(r) with periodicity T of h(t) so that just one sharp peak appears at 
h = —h in p{h(r)). This case would yield a rectangular hysteresis loop and hence with 
the largest possible area. On the other hand, if the passages take place all over randomly, 
for the other extreme case of least synchronization, p(h(r)) would be uniform resulting in 
a hysteresis loop of zero area. The hysteresis loop area, thus, provides a measure of the 
degree of synchronization of passages. It is to be noted that the calculation of hysteresis 
loop automatically takes into account of informations buried in all the peaks of the passage- 
time distribution p{r). It has earlier been shown H that hysteresis loop area so calculated 
acquires a maximum as a function of the input noise strength D as well as a function of the 
sweep rate \h\ (equivalent to the frequency) of the saw-tooth type periodic input signal. We 



now discuss our numerical experiment wherein the above explanation about the measure of 
synchronization should carry through, though not in as obvious a manner, from the residence 
time distributions in each of the two wells. 

In order to obtain the residence time distributions, pi(r) and p2(j) in each of the wells 
1 and 2, respectively, we monitor the trajectory m(t) of the particle [Eqs. (1) and (2)] for a 
long time and keep putting markers on the time axis whenever a passage from one well to 
the other takes place. We take h(t) of the same form as in Fig. Id but with h(t = 0) = 0. 
We consider the passage to take place only if the trajectory m(t) crosses the inflexion point 
on the far side of the maximum of the potential barrier separating the two wells. From the 
markers on the time axis we also obtain the jump field values h(t) for switching from one 
well to the other and hence the corresponding distributions pi 2 (h(t)) and p 21 (h(t)). Figures 
2a and 2b show the typical plots of pi(r) and P2i(h), respectively. Notice that we have, 
now, calculated pi 2 (h), etc., properly so that the distributions spread over the entire period 
of h(t), and includes both ascending as well as descending parts of h(t). It is now easy to 
calculate the probability m 2 (h) that the trajectory lies in the well 2 when the field value 
h(t) = h, from the discrete equation, 

m 2 (h) = m 2 (h - Ah) - m 2 (h - Ah)p 21 (h)Ah + m^h - Ah)p 12 (h)Ah, (4) 

with negligibly small Ah, and similarly for mi(h) = 1 — m 2 (h). We, then, calculate the 
hysteresis loop m(h) = m 2 {h) — mi(h) iteratively that satisfies the closed loop condition 
m(h(t + To)) = m(h(t)). Figure 2c shows a typical (asymptotically stationary) hysteresis 
loop. Please note that the hysteresis loop need not be saturated for every case in contrast 
to Fig. Id. We, now, discuss the results on the stochastic synchronization(SS) of passages 
from the variation of hysteresis loop area as a function of noise strength D and also as a 
function of sweep rate \h\. 

Throughout our calculation, we take a = 2.0 and b = 1.0 in the expression for U(m), so 
that the barrier height when h(t) = is 1. We take ho = .7h c and .9h c and calculate for 
each of these cases, the distributions pi(t) and p 2 (t) and also the passage field distributions 
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Pi2 anc ^ p2i(h(t)) for various values of h. The results obtained for both the field 
amplitudes 0.7h c and 0.9h c have qualitatively similar trends. Figure 3 shows the variation 
of hysteresis loop function of noise strength D. We observe that the hysteresis 

loop area initially increases as D is increased from a small value, attains a maximum value 
at D — D ma , say, and then decreases gradually as D is increased further. It shows that at 
D = D ma the output signal is most synchronized with the input signal h(t). The hysteresis 
loop area is also found to show maximum as a function of h (Fig. 3b). 

From the residence time distributions we calculate the Fourier transform and each of 
the components are squared to obtain the power spectral density ||. As expected large 
peaks are obtained at regular intervals. We calculate the ratio of the height of the first peak 
and the (background) noise level at the same frequency. This signal-to-noise ratio (SNR) is 
plotted in Fig. 4 as a function of D. The errors in calculating the SNR are large mostly 
because of the arbitrariness in fixing the (low) background noise level. However, the errors 
do not affect the general trends of our results. We, indeed, find SR in the usual sense of the 
SNR maxima. Interestingly, the value of D — DsNn{h) at which SR occurs, is not the same 
as the value of D — D ma (h) at which the passages are the most synchronized. However, 
DsNR^h) and D ma (h) are quite close and they tend to become closer as h — > 0, as shown in 
Fig. 5. (A finer examination to smaller h is beyond the computing power available to us.) 
SR and SS are, thus, not unrelated [0,H- SR may ,therefore, be a natural manifestation of 
SS. Moreover, SR, defined as the maximum of SNR appears not only as a function of noise 
strength D but also as a function of the sweep rate h (or equivalently the frequency) of the 
external field (input signal) as shown in Fig. 6 (cf. with Fig. 3b). 

In conclusion we state that hysteresis loop area, which is an average effect of all the peaks 
in the residence time distribution, is a good measure of synchronization of passages with the 
input signal. In an earlier work 0, by taking the only first few peaks, and in particular the 
first peak, of the residence time distribution, SR was given an alternative description and 
was shown to be a bonafide resonance. On closer examination we find that rest of the peaks 
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too play important role. We, however, concretise their assertion, by taking into account all 
the peaks of the residence time distribution through the hysteresis loop area, that SR is a 
genuine resonance, at least in a double-well potential system, and is due to the synchronized 
response of the system to the input periodic signal. Also, the degree of synchronization of 
passages and the SNR show maximum as a function of the sweep rate (or equivalently the 
frequency) of the input signal |§. 
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FIG. 1. Shows (a) the time variation of the external field h(t), (b) the first-passage-time dis- 
tribution p(r), (c) the passage field distribution p(h), and (d) the corresponding hysteresis loop 
calculated from p(h), for ho = 0.7h c , D = 0.3, and the period of h(t), Tq = 28.0. 

FIG. 2. Shows (a) the residence time distribution pi(r) in the well 1, (b) the passage field 
distribution p2i(h) for passage from well 2 to well 1, and (c) the corresponding hysteresis loop 
m(h) for h = 0-9h c , D = 0.2, and period of h(t), T = 36.0. 

FIG. 3. Plots of (a) hysteresis loop area A as a function of D for h = 0.05/i c (o), 0.1/i c (D), 
0.2/i c (o), 0Ah c (A), and 0.6/i c (v) 5 and (b) hysteresis loop area A as a function of field sweep rate 
h for D = 0.1(o), 0.15(a), 0.2(o), and 0.3(A) for h = 0.9h c . In this and in the rest of the figures 
the lines joining the points are only to guide the eye. 

FIG. 4. Plots of signal-to-noise ratio (SNR) as a function of D, for h/h c = 0.05(o), 0.1(D), 
0.2(o), 0.4(A), 0.5(<l), 0.6(v), and 0.72(>), for ho = 0.9h c . 

FIG. 5. Plots the peak positions of the plots of hysteresis loop area A versus D(o) and those 
of SNR versus D(D), for ho = 0.7h c (empty symbols, solid joining lines) and for ho = 0.9h c (filled 
symbols, dashed joining lines). 

FIG. 6. Shows SNR as a function of h for ho = 0.9h c , for various values of D = 0.1(o), 0.15(D), 
and 0.2(o). 
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